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A high level polarizable force field is used to study the temperature dependence of hydrophobic 
hydration of small-sized molecules from computer simulations. Molecular dynamics simulations 
rS I of liquid water at various temperatures form the basis of free energy perturbation calculations 
I that consider the onset and growth of a repulsive sphere. This repulsive sphere acts as a model 
construct for the hydrophobic species. In the present study an extension is pursued to all in 
all seven independent target temperatures starting close to the freezing point and ranging up 
to almost the boiling point of liquid water at standard conditions. Care is taken to maintain 
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proper physico-chemical model description by cross-checking to experimental water densities 
at the selected target temperatures. The polarizable force field description of molecular water 
turns out to be suitable throughout the entire temperature domain considered. Derivatives 
of the computed free energies of hydrophobic hydration with respect to the temperature give 
access to the changes in entropy. In practice the entropy differential is determined from the 
negative of the slope of tangential lines formed at a certain target temperature in the free 
energy profile. The obtained changes in entropy are of negative sign for small sized cavities, 
hence reconfirm basic ideas of the Lum Chandler Weeks theory on hydrophobic hydration of 
small-sized solutes. 

1 Introduction 

The hydrophobic effect is widely believed to play a decisive role in protein folding, one of the 
key challenges in biophysical science and research of today lU [21 [3]. Theoretical studies on 
the hydrophobic effect are of great relevance to a broad range of biosciences and a deeper 
understanding of hydrophobicity could certainly have its beneficial infiuence on many central 
questions in current biophysical research. Among others, Lum Chandler Weeks (LCW) the- 
ory of hydrophobicity [1] has received widespread appreciation. LCW theory describes the 
hydrophobic effect in terms of reorganizational work due to maintainance of a hydrogen bond 
network established between individual water molecules. A difference is made between small- 
sized hydrophobic solutes (volume dependence) and large-sized hydrophobic solutes (surface 
area dependence). Water molecules are believed to re-arrange appropriately around small-sized 
hydrophobic solutes thereby inducing the formation of clathrate-like substructures without 
destruction of the hydrogen bond network. In contrast, large-sized hydrophobic solutes are 
thought to enforce a complete re- arrangement of the hydrogen bond network adjacent to the 
hydrophobic solute, which in turn leads to agglomeration, aggregation and precipitation of the 
large-sized hydrophobic molecules. Such a picture of hydrophobicity would render the change 
in entropy for the process of solvating small-sized hydrophobic solutes to be negative in sign 
(increase in order), and the opposite for large-sized hydrophobic molecules (decrease in order). 
Hence, two immediate questions arise naturally: i) can computer simulations verify the claim of 
entropic drop for hydration of small-sized hydrophobic solutes ? ii) are current model descrip- 
tions of molecular water able to reproduce the physics of hydrophobicity correctly ? Although 
complementary, these two questions have to be addressed by atomistic computer simulations 
of hydrophobic solutes in aqueous solution. 

Computer simulations have become a valuable tool in the study of hydrophobicity O O [3, [3 [9] . 
Among other techniques, free energy calculations have been introduced and advanced to directly 
study AG trends for various physico-chemical processes [IHl [HI • One such possible process 
is to investigate the AG corresponding to the introduction, onset and growth of a repulsive 
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sphere located in the center of a simulation cell filled with water molecules. In such a model, the 
repulsive sphere stands as a representative of an artificial hydrophobic solute. The associated 
free energy change is known as the cavitation free energy, AG™'' [13]. Postma and coworkers 
have introduced this type of calculation as one of the early examples of Free Energy Perturbation 
calculations (FEP). Their approach is known as the Overlapping Spheres Technique (OST) [TT] . 
Recent re-evaluations and variations of the OST in the context of hydrophobicity have been 
summarized in [T3] . The advantage of true estimates of AG becomes evident when looking at 
the derivative with respect to temperature, ^(AG) = ^{AH — TAS) = —AS, i.e. from a 
record of different AG values at different temperatures the change in entropy can be determined 
from the negative slope of the tangent formed at a particular point. Thus if the LOW arguments 
hold true, then for solvating small-sized hydrophobic solutes one would obtain a bell-shaped 
curve of the temperature plot of AG, for otherwise the slope of the tangent at room temperature 
can not become positive, hence AS not negative. However, high temperature simulations in the 
molecular mechanics description need to be carefully cross-checked against experimental data 
whether they still can provide reasonable representations of the native state. Therefore, only 
if one can confirm that a set of high temperature simulations still takes place at reasonable 
physico-chemical conditions, an estimate of entropies can be obtained. For the type of atomistic 
interaction governing hydrophobic phenomena, a proper account of the liquid water density at 
high temperatures is the most essential precondition for pursuing estimates of the entropy 
[I51[T6]. 

Model descriptions of water have a great infiuence on the outcome of biomolecular simulations 
p!7l [T8] but the employment of prominent water models EQ] is common practice in present 
biophysical research. A specific high level description of molecular water has been proposed 
with the use of polarizable models [21]. The just cited AMOEBA model has been shown 
to exhibit excellent description of the temperature and pressure dependence of water [22] . 
AMOEBA water was successfully applied in describing aspects of the hydrophobic effect [23] . 
Since this previous study was focussing on room temperature behavior, an extension towards 
high temperature repeats is straightforward as it could lead to new insight into fundamental 
principles of hydrophobicity. Other studies on the entropy of hydrophobic hydration were 
mainly based on rigid water models so far [2H [25] . 

The present article reports cavitation free energy calculations performed at seven individual 
temperatures in the range of 277 K to 370 K. The polarizable AMOEBA water model [21] 
is used and Ewald summation is applied within the Molecular Dynamics simulations (MD) 
that form the basis of the FEP calculation of the OST approach. Emphasis is placed on 
verification of proper physico-chemical model description at high temperatures by comparison 
to experimentally obtained trends of liquid water densities. Entropic changes for this process 
are derived from the temperature dependence, that is the entropy of hydrophobic hydration is 
estimated by means of computer simulations. The computational demand of this study is on 
the order of three CPU years on decent architectures (Itanium 1.4 GHz) [26] and can only be 
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satisfied from massive employment of grid computing systems, such as for example the Austrian 
grid [27]. 
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Methods 



2.1 Simulation Cell Set Up 



Seven individual cubic boxes composed of 6 x 6 x 6 grid cells were formed, where the sub- volume 
of the grid cells was adjusted to reproduce the experimental liquid water density corresponding 
to a chosen target temperature. Seven target temperatures were selected, 277 K, 300 K, 315 
K, 330 K, 345 K, 365 K and 370 K. The structure of a single water molecule was optimized 
and then periodically translated and copied to each of the centers of the grid cells. Thus all 
simulation cells contained 216 water molecules. After initial construction the systems were 
minimized and subjected to simulated annealing using 2000 steps of 1.0 fs each to approach 
1000 K peak temperature before linearly cooling down to one of the seven target temperatures. 
Volume modification during simulated annealing was less than 1 % when compared to the 
box dimensions upon start-up. All calculations were performed with the TINKER package for 
molecular modeling version 4.2 [33] . 




24 individual MD-trajectories per chosen target temperature (see section 12. ip were recorded. 
The TINKER package [33| was used for all computations. Polarizable AMOEBA force field pa- 
rameters for water [21] were employed. AMOEBA works on the basis of self-consistent induced 
atomic dipoles. NpT ensembles were selected at 1.0 atm target pressure with the tempera- 
ture/pressure coupling method due to Berendsen [31] which accounts for box-size changes when 
the volume of the repulsive cavity is introduced. Default coupling constants to the thermostat 
and barostat were used, that is 0.1 ps and 2.0 ps respectively. Other choices could have been 
made [35] , especially with respect to the "Flying ice cube" phenomenon [36] , but here we tried 
to follow closely the protocol used during AMOEBA development. Time steps of 1.0 fs were 
employed, no restraints/constraints were applied and Ewald summation was used. Two types 
of perturbation potentials were introduced. 




2.2 MD/FEP Calculations 
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to avoid discontinuities when the thermal radius of the repulsive cavity, B*, approaches A 
[371 [38] • Parameter A in equations ([1]) and ([2[) describes the degree of perturbation and assumes 
values from A = (unperturbed) to A = 1 (fully perturbed). All technical details concerning 
the perturbation potential are given in [28]. Trajectories were recorded for 100 ps each and the 
perturbations shown in equations ([T]) and ([2]) were calculated and stored every time step. As 
reported previously, the A = 0.5 simulation applying V^^^'^ of equation ([2]) is quasi-unperturbed 
and was used to monitor the time evolution of the macroscopic liquid water density. Data 
evaluation was started only after stable levels of the liquid water density had been reached 
(approximately after 40 to 50 ps, see Tabled]). The A = 0.5 simulation has also been analyzed 
with respect to maintainance and changes in the Hydrogen Bond network, as outlined in [H] |12] 
(see Supplementary Material). 

2.3 Data Evaluation 

For each of the 7 selected target temperatures (see section 12. ip the recorded 25 trajectories 
were used for application of Zwanzig's [10] formula 

AGi\i) = -kBT /n^e"^f^^^'+'^^^'^^^'^^^ (3) 

with < > depicting a thermodynamic average, ks being the Boltzmann constant and Ti. the total 
energy of the system. The average is formed at Aj with small perturbations 6X around Aj. The 
OST is applied similarly to the description given in [28]. Care has been taken of using proper 
ksT factors for all the seven different temperatures mentioned above. This will effect equation 
([3]) and the conversion from repulsive radii, B, to thermal radii, B*, i.e. B* = B{lkBT)~T2 . 
Perturbations exceeding energies of 2k bT were not considered for total averages (PEP require- 
ment). There are however always sufficiently enough alternative overlap combinations in all 
the simulations to smoothly connect two adjacent repulsive spheres and get statistical averages. 
Obtained raw data from the OST were fitted with polynomials of degree 2 and corresponding 
coefficients k^, ki and are summarized in Table [H Resulting cavitation free energies from 
the coefficients in Table [T] are expressed in units of kcal/mol (1 kcal = 4.184 kJ). 

3 Results 

3.1 AMOEBA based computer simulations of molecular liquid wa- 
ter confirm LCW ideas of the unit volume entropic change for 
solvating small-sized hydrophobic solutes 

MD/FEP calculations based on the polarizable force field AMOEBA for the computation of 
cavitation free energies in liquid molecular water are carried out at seven different target tem- 
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peratures. The results are fitted similarly to previous calculations [2S] and resulting coefficients 
ko, ki, k2 are summarized in Table [1] {revised Pierotti Approach, rPA). The coefficients are 
used to calculate temperature trends for growing cavities of perfectly spherical shape. Figure 
[1] shows a comparison of these data in reduced energy units normalized to unit volumes. This 
type of reduced units allows one to immediately compare unit volumes of 1 cubic A to each 
other. Unit volume cavitation free energies of small-sized cavities (squares and discs in Figure 
[1]) clearly exhibit a positive slope for the tangent formed at room temperature data points. This 
indicates a drop in entropy for the process of hydrating a hydrophobic volume of 1 cubic A. 
Similar constructions of tangents for larger-sized cavities (triangles and diamonds in Figure [1]) 
show that the positive slopes become smaller as cavities grow. Thus high level force field based 
MD/FEP calculations verify the anticipated change in unit volume entropy following LCW 
theory for small-sized hydrophobic solutes (see for example Figure 2 in [29] for comparison to 
Figured] of the present work). 

3.2 Cavitation entropies of small sized hydrophobic solutes decrease 
steadily with increasing cavity size 

Equally interesting than unit volume entropies are the changes observed for full size cavitation, 
that is trends due to the creation of the entire cavity volume. Figure [2] shows the temperature 
trend of cavitation free energies of small sized hydrophobic solutes without scaling to the 
unit volume. Tangential lines formed at room temperature data points again exhibit positive 
slopes thus again indicating a drop in entropy. However, the magnitude of the slope becomes 
larger with increasing cavity size (e.g. compare slopes of tangential lines for small-sized cavities, 
squares and discs in Figure[21 to tangential slopes of larger-sized cavities, triangles and diamonds 
in Figure [2]). Inversion in the sign of the slopes of tangential lines is seen close to the boiling 
point. Quantification of individual enthalpic and entropic contributions is given in Table [2l 
Small negatively signed values of AH are likely to result from numerical processing and should 
rather not taken to represent physical meaning. As may be seen from Table [2], the process is 
largely dominated from entropy and only to a minor extent enthalpic. Therefore we need to 
conclude that the change in entropy is steadily decreasing the larger the hydrophobic solute 
becomes. 

3.3 Elevated temperature simulations of molecular liquid water 
largely resemble the experimentally observed temperature de- 
pendence of the liquid water density 

In order to assure proper model description of intermolecular forces at temperatures other 
than room temperature, the macroscopic liquid water density is extracted from each simulation 
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performed at a certain target temperature. Figure El shows a comparison between simulation 
data (squares) and experimental measurements [30] (discs). Largest deviations occur at 300 K 
and 370 K with a maximum unsigned error of 0.01 appearing at 300 K. The root mean 
square deviation amounts to 0.0056 Given the rather close match at temperatures close 
to the boiling point, the overall rating of the simulations concerning intermolecular interac- 
tions at elevated temperatures must be considered very satisfactory. Consequently, taking the 
present data to form T-derivatives in order to determine entropic changes seems to be a valid 
approximation. 

3.4 Rather close quantitative agreement between LCW predictions 
of the unit volume entropic change for the solvation of small- 
sized hydrophobic spheres and results derived from present com- 
puter simulations 

Graphical extrapolation of unit volume AS values from the slopes of tangents formed at room 
temperature data points in temperature dependence plots of the unit volume AG'^'**' (Figure [1]) 
yields the following values: AS » ^ -0.00027 ^E^^B^, AS » ^ -0.00020 iE^A^Hi , ^ 

B=2A ^ ' B=3A ^ ' B=4A 

-0.00015 '^'^t'a™"^ and AS ° ~ -0.00012 i^£^/i££l. \ similar evaluation of the data presented in 

A K B=5A ^ 

l29l leads to comparable values, i.e. aS^"'^"?''"'^'^'^'^'°'^ ~ -0.00007 '^'^t'a . Thus present computer 

B=3A A K 

simulation data not only show qualitative agreement with LCW theory on the hydrophobic 
effect of small-sized solutes, but also lead to comparable results in absolute numbers of unit 
volume entropic changes. 

3.5 Present computer simulations are limited to the small-size do- 
main of hydrophobic hydration 

rPA coefficients obtained from MD/FEP calculations in AMOEBA water are strictly valid only 
in the domain of existing FEP data. It is reasonable to extrapolate into the extended cavity size 
domain due to the smoothness of the data. However, after about B=5 A (twice the radius of 
the largest accumulated perturbation) any further usage is certainly speculative. Therefore it 
was surprising to see indications of general applicability even in largely extended cavity domains 
|14j . In order to probe the quality of large scale extrapolation the temperature trend of the 
surface tension, a, is shown in Figure [Hand compared to experimental data [31]. Approximation 
of a from the present rPA data is critical (see Discussion). It involves limit value consideration 
of lim A G™''/(47r B^) [Tl] and the presented data is due to setting B=100 A, hence clearly 

B— >oo 

an extrapolation to large scale. Figure H] shows that the calculated surface tensions are close to 
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the experimental values but do not reproduce the correct trend with increasing temperature, 
demonstrating the anticipated uncertainty with large-size approximations. 

3.6 Analysis of the hydrogen bond network as determined by the 
ensemble of snapshot structures obtained from AMOEBA based 
MD simulations goes hand in hand with recently reported find- 
ings of experimental as well as theoretical studies 

MD/FEP runs carried out at A = 0.5 (quasi-unperturbed) are examined with respect to struc- 
tural relationships maintained between individual water molecules. Corresponding sets of saved 
snapshot structures of individual water boxes at T=300 K and T=370 K are analyzed. The 
three basic Radial Distribution Functions (RDF) that characterize the "fine structure" of bulk 
water are derived and corresponding trends shown in the Supplementary Material, Figures A-C 
(300 K) and Figures D-F (370 K). RDFs from the present data reproduce those reported in the 
original AMOEBA paper (see [21] Figures 5-7) very well. The RDF of the oxygen-oxygen 
distance, goo(r), shows the characteristic immediate neighborhood peak at 2.9 A and a second 
very diffuse peak around 4.6 A which is hardly distinguishable from the baseline. The second 
peak in the goH(i') — due to pairs of water molecules directly associated via hydrogen bonds 
(HBs) — appears at 1.9 A in the present data. The gHH(r) exhibits a characteristic second 
peak at 2.5 A. Following the procedure devised in [12] for geometrical constraints defining a 
particular HE we determine an average number of established HBs of 1.6 per water molecule 
in the 300 K simulation data. The frequency of tetrahedral coordination mediated by HBs is 
only on the order of 2 % at this simulation temperature. This seems to be consistent with 
the experimental study of Wernet et al jH] where the conclusion was drawn that bulk water 
predominantly exists in a state with 2 HBs established per water molecule. The present data 
is also in close agreement with the MP2 results of the QM/MM study by Xenides et al [42] . 
Significant changes observed at the high temperature simulation of T=370 K include the entire 
loss of the second broad shoulder around 4.6 A in the goo(i')? the decrease of the average 
number of HBs to 1.4 per water molecule and the general loss of tetrahedral coordinations 
mediated by HBs. 

4 Discussion 

The present study has employed the high level polarizable AMOEBA force field for computer 
simulation of liquid water [21] at several temperatures and normal pressure. We are interested 
in the free energy cost of creation of small-sized cavities, which are model constructs for hy- 
drophobic molecules. A first requirement was to show that intermolecular relationships are still 
maintained at reasonable physical conditions upon temperature increase A rather direct 
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evidence in this regard is the relatively close match between simulated and experimental data of 
the liquid water density shown in Figure [31 Next the temperature profile of AG was used to get 
estimates for the unit volume AS (Figured]) as well as for the full size AS of cavitation (Figure 
[2]). Computer simulation values for AS can be derived from the negative slope of the tangents 
constructed at a certain temperature. In so doing the AS for room temperature hydrophobic 
hydration of small-sized molecules was shown to exhibit a negative sign, which was predicted 
previously from LCW-theory [3]. The absolute value of the change in cavitation entropy, AS, 
becomes larger with increasing size of the hydrophobic molecule (see Figure |2] and Table [2]). 
Tendency inversion is observed only at elevated temperatures close to the boiling point. 
In general, the LCW picture of the unit volume hydrophobic effect of small-sized solutes is 
very well reproduced (see for example Figure 2 in [29]). Small differences include, the exact 
position of the maximum, which occurs closer to the boiling point in the present data, the 
quantification of the unit volume AS, which leads to higher values in the present study and 
the continuity of the onset of the curves at temperatures close to the freezing point, which is 
steady here, but shows some inversion at B=4 A solutes in LCW-theory. Extrapolation of the 
present data to very large solute sizes is impossible as seen from the wrong prediction of the 
temperature profile of the surface tension (Figure H]). The present approach has made use of 
a very expensive computational technique (Ewald sum MD/FEP within the OST) and could 
only be carried out under massive employment of grid computing [27|. All final rPA coefficients 
corresponding to all the considered temperatures are summarized in Tabled] It is interesting to 
note, that all present kQ, ki, /c2 coefficients derived for 300 K lead to only minor alterations in 
the data presented recently in [T3] that were based on non-Ewald simulations, but application 
of plain periodic boundary conditions instead [23] . 

An interesting observation made in the present work is the fact that the overall entropy of 
hydrophobic hydration is steadily decreasing with growing solute size (see Figure [2] and Table 
[2]). This is somewhat contradictory with the LCW-picture of the anticipated change in the 
sign of AS when the solutes are thought to cross over from small length scales (volume de- 
pendence) to large length scales (surface area dependence). The critical crossover dimension 
was approximated to be on the order of 10 A of radial extension of the hydrophobic solute [1]. 
On the other hand, strictly speaking, the crossover region is beyond the reach of the current 
data set and firm statements can only be made upon extension of the present approach into 
the critical domain of crossover length scales. This would however involve an even bigger ini- 
tiative of super computing with unforeseeable complexity. It is still interesting to note that an 
extension of the rPA approach into medium length scales was recently shown to be not entirely 
unreasonable when operating with the solvent excluded volume |14j . 

Predicting the temperature dependence of the surface tension (see Figure H]) from present rPA 
data failed. First of all it should be noted, that surface tension per se is appropriate for the 
description of water embedded macroscopic bubbles. These bubbles are fundamentally different 
to cavities of comparable size, because the latter are strictly empty. Therefore, all what makes a 
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bubble stabilize, i.e. the vapour-like molecules in the interior bouncing back and fourth against 
the bubble walls, will be completely missing in a cavity, because the cavity interior is void. 
Nevertheless, the equivalence in units of surface-normalized cavitation free energies on the one 
side and surface tensions on the other side, makes it attractive to try the prediction of the 
latter via lim A G'^"^/(47r B^). However, the energy in the context of surface tension implies 

an energy required to maintain the liquid-vapour coexistence of a bubble, while it expresses an 
amount of clearance work to create the empty space in a cavity. These two different types of 
energies should not be confused with each other and in principle, need not even be correlated. 
Thus the rather successful match of lim A C^""/ (47r B^) with the experimental surface tension 

B^oo 

of water reported in [H] could also be purely fortuitous. On the other hand it could also explain 
why surface tension based solvation models [32] and LCW-theory have become so successful. 
LOW theory — like the earlier concepts introduced by Kauzmann [39] and Tanford [ID] — 
build upon oil/water surface tension and explain the switch in sign of AS" for the hydration of 
hydrophobic solutes of increasing size by the experimental temperature dependence of surface 
tensions. LOW in particular has proclaimed that such a switch in sign would occur at spherical 
volumes of radial dimension on the order of 10 A. Taking up this idea, AS* values in the vicinity 
of the 10 A domain should smoothly become smaller and smaller before they reach zero and 
then start to assume negative values. A plot of AS* in this critical region should therefore 
certainly not exhibit characteristics of a monotonic growth, otherwise the transition would 
have to occur discontinuously. In the present study, although still operating far away from the 
critical domain of 10 A, we do however observe signatures of a monotonic growth (see Table 
[2]). This will be of particular interest to examine in closer detail from extended studies of the 
type presented here. 



5 Conclusion 

In summary, evidence has been presented that the hydrophobic effect of small-sized solutes is 
due to the drop in entropy. The insight gained here is also an encouraging sign that high level 
force fields are approaching now standards that allow one to study effects beyond the reach of 
experimental methods. 
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Table 1: Summary of the data resulting from free energy perturbation calculations of repul- 
sive spheres in water at various simulation temperatures. Equilibration values of liquid water 
densities (second column) are derived from A = 0.5 runs (quasi- unperturbed) after neglecting 
an initial fraction of the data set (fourth column). The data are fitted with polynomials of 2'^'^ 
degree and resulting rPA-coefficients ko, ki and are included in columns 5-7. 



5cxp Initial rPA-Coefficients 



Discard ko ki /c2 

il / mol 1 r kcal / 1 



[K] [i] [;^] [ps] [kcal/mol] [^] [^^1 



277 


0.995 


1.000 


40 


0.456 ± 0.015 


-1.602 ± 0.062 


1.116 ± 0.082 


300 


1.006 


0.996 


40 


0.427 ± 0.012 


-1.594 ± 0.044 


1.183 ± 0.046 


315 


0.997 


0.992 


40 


0.418 ± 0.041 


-1.613 ± 0.060 


1.220 ± 0.022 


330 


0.988 


0.985 


40 


0.398 ± 0.008 


-1.607 ± 0.045 


1.236 ± 0.052 


345 


0.980 


0.977 


40 


0.405 ± 0.011 


-1.607 ± 0.045 


1.251 ± 0.045 


365 


0.963 


0.964 


40 


0.383 ± 0.015 


-1.580 ± 0.026 


1.262 ± 0.038 


370 


0.953 


0.960 


50 


0.411 ± 0.009 


-1.611 ± 0.037 


1.236 ± 0.029 
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Table 2: Decomposition of cavitation free energies AG'^'*'^ (column to the right) of small-sized 
spherical volumes of radius B (leftmost column) into entropic (third column) and enthalpic (sec- 
ond column) fractions. The defragment at ion is based on initial finite difference determination 
of AS from the temperature trend of AG'^"'^ (see Figure [2]) with subsequent trivial derivation 
of AH. The analysis is made for room temperature behaviour (T=300 K). 



B 


AH 


-TAS 




[A] 


[kcal / mol] 


[kcal / mol] 


[kcal / mol] 


2 


-0.728 


2.699 


1.972 


3 


-0.335 


6.628 


6.293 


4 


0.798 


12.183 


12.981 


5 


2.672 


19.364 


22.036 
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1 Temperature dependence of the unit volume free energy of cavitation for small 
repulsive spheres of growing size. Small-sized cavities of radii B = 2 A (squares) 
and B = 3 A (discs) exhibit room temperature tangents with positive slope that 
correspond to a unit volume AS of negative sign. Similar analysis of tangential 
slopes for larger volumes, i.e. B = 4 A (triangles) and B = 5 A (diamonds) 
demonstrates that the magnitude of the slopes decreases at 300 K when the size 
of the solute becomes bigger 

2 Temperature dependence of the full size free energy of cavitation for small repul- 
sive spheres of growing size. Room temperature tangents exhibit positive slopes 
that correspond to a AS of cavitation of the entire cavity volume with negative 
sign. The magnitude of the slope of tangential lines increases with growing cavity 
size, i.e. when following cavities of radii B = 2 A (squares), B = 3 A (discs), B 
= 4 A (triangles) and B = 5 A (diamonds). Inversion in the sign of the slopes 
of tangential lines is seen close to the boiling point 

3 Comparison of simulated and measured macroscopic density, p, of liquid water at 
various temperatures. The simulated water densities (squares) are derived from 
A=0.5 MD/FEP calculations that take place quasi-unperturbed [23]. The exper- 
imental densities (discs) are obtained from standard tabulations [3D]. Greatest 
deviations from experimental values occur at 300 K and 370 K. The root mean 
square deviation is 0.0056 gcm~^ 

4 Comparison of simulated and measured surface tension, a, of liquid water at 
various temperatures. The calculated values of a (squares) are due to lim A 

B— >oo 

G™''/(47r B^) with B=100 A. Experimental reference data (discs) have been 
obtained from standard tabulations [31]. Although qualitatively comparable, 
the calculated values of a exhibit inverse temperature profiles, which might be 
an artefact stemming from large-scale extrapolation 
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Figure 1: Temperature dependence of the unit volume free energy of cavitation for small 
repulsive spheres of growing size. Small-sized cavities of radii B = 2 A (squares) and B = 
3 A (discs) exhibit room temperature tangents with positive slope that correspond to a unit 
volume AS of negative sign. Similar analysis of tangential slopes for larger volumes, i.e. B 
= 4 A (triangles) and B = 5 A (diamonds) demonstrates that the magnitude of the slopes 
decreases at 300 K when the size of the solute becomes bigger. 
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Figure 2: Temperature dependence of the full size free energy of cavitation for small repulsive 
spheres of growing size. Room temperature tangents exhibit positive slopes that correspond to 
a AS of cavitation of the entire cavity volume with negative sign. The magnitude of the slope 
of tangential lines increases with growing cavity size, i.e. when following cavities of radii B = 
2 A (squares), B = 3 A (discs), B = 4 A (triangles) and B = 5 A (diamonds). Inversion in 
the sign of the slopes of tangential lines is seen close to the boiling point. 
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Figure 3: Comparison of simulated and measured macroscopic density, p, of liquid water at var- 
ious temperatures. The simulated water densities (squares) are derived from A=0.5 MD/FEP 
calculations that take place quasi- unperturbed [23]. The experimental densities (discs) are ob- 
tained from standard tabulations [30]. Greatest deviations from experimental values occur at 
300 K and 370 K. The root mean square deviation is 0.0056 gcm~^. 
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Figure 4: Comparison of simulated and measured surface tension, a, of liquid water at various 
temperatures. The calculated values of a (squares) are due to lim A G'^™/(47r B^) with 

B— »oo 

B=100 A. Experimental reference data (discs) have been obtained from standard tabulations 
|31j . Although qualitatively comparable, the calculated values of a exhibit inverse temperature 
profiles, which might be an artefact stemming from large-scale extrapolation. 
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